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ABSTRACT 


Context. Polarization of radio emission in extragalactic jets at a sub-milliarcsecond angular resolution holds important clues for 
understanding the structure of the magnetic field in the inner regions of the jets and in close vicinity of the supermassive black holes 
in the centers of active galaxies. 

Aims. Space VLBI observations provide a unique tool for polarimetric imaging at a sub-milliarcsecond angular resolution and 
studying the properties of magnetic field in active galactic nuclei on scales of less than lO"* gravitational radii. 

Methods. A space VLBI observation of high-redshift quasar TXS 0642-1-449 (OH 471), made at a wavelength of 18 cm (frequency of 
1.6 GHz) as part of the Early Science Programme (ESP) of the RadioAstron mission, is used here to test the polarimetric performance 
of the orbiting Space Radio Telescope (SRT) employed by the mission, to establish a methodology for making full Stokes polarimetry 
with space VLBI at 1.6 GHz, and to study the polarized emission in the target object on sub-milliarcsecond scales. 

Results. Polarization leakage of the SRT at 18 cm is found to be within 9% in amplitude, demonstrating the feasibility of high fidelity 
polarization imaging with RadioAstron at this wavelength. A polarimetric image of 0642-t449 with a resolution of 0.8 mas (signifying 
an ~ 4 times improvement over ground VLBI observations at the same wavelength) is obtained. The image shows a compact core¬ 
jet structure with low (»! 2%) polarization and predominantly transverse magnetic field in the nuclear region. The VLBI data also 
uncover a complex structure of the nuclear region, with two prominent features possibly corresponding to the jet base and a strong 
recollimation shock. The maximum brightness temperature at the jet base can be as high as 4 x lO'^ K. 


1. Introduction 

Very long baseline interferometry (VLBI) observations in which 
one of the antennas is placed on board of an Earth satellite (space 
VLBI) have provided a capability of reaching unprecedentedly 
high angular resolution of astronomical observations ( |Arsentev| 
et al.|1982| Levy et al.|1986[ Hirabayashi et al.|2000| l. The space 
VLBI mission RadioAstron combines ground-based radio anten¬ 
nas operating at frequencies 0.32, 1.6, 5, and 22 Gigahertz (GHz) 
with a 10-meter antenna (Space Radio Telescope, SRT) on board 


of the satellite Spektr-R launched on 18 July 2011 (Kardashev 
et al.|2013k 

At 0.32, 1.6, and 22 GHz, the SRT delivers data in both left 
(LCP) and right (RCP) circular polarization at 0.32, 1.6, and 
22 GHz, which enables fully reconstructing the linearly polar¬ 
ized emission from target objects at these frequencies. At 5 GHz, 
the SRT provides only LCP data because of a failure of the on¬ 
board hardware. 

After its launch, the SRT was first tested in a single-antenna 


mode during the in-orbit-checkout (IOC) period (Kovalev et al. 
|2014| l and later in the interferometry mode, in combination with 
ground-based antenna (fringe search, PS), successfully deliver¬ 


ing fringes at each of the four observing bands (Kardashev et al. 

[2fe][2nT4al l. 

The ensuing early science program (ESP) of RadioAstron 


started in Pebruary 2012 and continued through June 2013 (Kar- 


Idashev et al.|2014a|l. One of the prime objectives of the program 


was to provide a bridge between the initial experimental mode 
of operations, observing, and data processing, and routine oper¬ 
ations that started after completion of the ESP. 

One of the main goals of the ESP was to establish the feasi¬ 
bility of imaging experiments with RadioAstron. The first imag¬ 
ing experiment with RadioAstron, targeting the bright and com¬ 
pact radio source 0716-1-714, was made during the ESP period at 
5 GHz (Kardashev et al. 2014a| l. Polarization imaging was tested 
with the observation presented in this paper. 

Polarization imaging at high angular resolution, both along 
and transverse to the jet direction, is crucial for measuring the 
magnetic field structures of jets to distinguish between differ¬ 
ent possible magnetic field configurations. The improvement of 
angular resolution provided by space VLBI observations makes 
it particularly interesting to compare the brightness temperature 
and polarization measurements made with RadioAstron to the 
ground VLBI with a similar resolution obtained at higher observ¬ 
ing frequencies. It has been shown that a minimum of several, in¬ 
dependent, resolution elements across a jet are necessary to pre¬ 
vent spurious results from being obtained ( |Hovatta et al.|201^ . 
Space VLBI observations are essential for resolving transversely 
a significant number of jets and studying the structure and prop¬ 
erties of their magnetic field. 

Leasibility of polarization measurements with space VLBI 
has been demonstrated with the VSOP ( Hirabayashi et al.|2000) l 
observations made with a satellite antenna recording data only in 
single (left circular) polarization (|Kemball et al.|2000{ [Gabuzda] 
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I& G6mez|2001| l. Results from the first polarization observation 
with RadioAstron at a frequency of 1.6 GHz are reported here. 
This observation was made in March 2013 as part of the ESP 
and it targeted a bright, compact radio source TXS 0642+449 
(OH471). 

The radio source 0642+449 is a distant low-polarization 
quasar (LPQ) located at a redshift of 3.396 ( |Osmer et al.|1994| l 
which corresponds to the luminosity distance of 30.2 Giga- 
parsecs (Gpc) and a linear scale of 7.58 kpc/arcsecond, assum¬ 
ing the standard ACDM cosmology ( [Planck Collaboration et al. 
|2015| l. The optical spectrum of 0642+449 does not show strong 
broad lines ( |Torrealba et al.|2012| l. In radio, its radio structure 
has a compact “core-jet” morphology (|Gurvits et al.||1992 Xu 


et al. 


1995|l, with a flat spectral index a ^ -0.1 ([Hovatta et al. 


2014 1 of the milliarcsecond (mas) scale core (bright narrow end 
of the jet). 

The maximum brightness temperature, Tb, measured in 
0642+449 decreases with observing frequency. Previous VSOP 
space VLBI observations o f 0642+449 at 5 GH z put a lower 
limit of Tb > 3.3 x 10^^K ( Dodson et al.||2008 l, while ground 
VLBI observations at 15 and 86 GHz resolved the compact core 
and yielded somewhat lower br ightness temperature s of 1.2- 
4.3 xl0‘^ K ([Kovalev et al.|2005| l and 1.1-1.6 xl0“ K ( |Lobanov 
et al.|2000l Lee et al.|2008|l. This implies that the compact core 


becomes gradually resolved out at higher frequencies, calling for 
comparison with the new space VLBI observations by RadioAs¬ 
tron. 

Ground VLBI monitoring of the source has revealed a range 
of proper motions of 0.012-0.078 mas/yr in the jet ([Lister et~ar 


2013| l and provided estimates of the jet Lorentz factor, Tj = 5.4, 


viewing angle 0, = 0.8° and intrinsic opening angle cpi - 0.3° 


(Pushkarev et al. [2009 1 . With these parameters, an equiparti- 
tion magnetic field of ^ 0.3 G can be estimated at the location 
of the jet core ( |Lobanov1|1998| l observed at 1.6 GHz, based on 
measurements of frequency dependent shift of the core position 
which have been reported for 0642+449 between frequencies of 
8.1, 8.4, 12.1, and 15.4GHz ( [Pushkarev et al.|2012D . 

Ground VLBI polarization observations at 5.0 and 8.4 GHz 
( [O’Sullivan et al.|2011[ l have shown a small (1.1-2.6%) degree 
of polarization in the core, with no conclusive evidence for Fara¬ 
day rotation (x - 15.6° and 15.3° at 5.0 and 8.4 GHz, respec¬ 
tively). At higher frequencies, between 8.4 and 15.1 GHz, a ro¬ 
tation measure (RM) of -280 rad/m^ has been reported in the 
core ( [Hovatta et al.|201^ . The Galactic RM in th e direction of 
0642+449 is estimated to be -5.4 + 18.3 rad/m^ (Taylor et al.' 
|2009[|0ppermann et al.[2012[ l. 

The polarization observations of 0642+449 with RadioAs¬ 
tron signal the first successful space VLBI polarization experi¬ 
ment using dual circular polarization measurements made at the 
space-home antenna. The methodology for RadioAstron polari- 
metric observations and data correlation is presented in Sect. ^ 
The RadioAstron observation of 0642+449 and reduction of the 
polarimetric data are described in Se ct. The polarization im¬ 
age of 0642+449 is discussed in Sect. |4| and compared with pre¬ 
vious ground VLBI observations of this object. 


2. Imaging with RadioAstron 

RadioAstron observations are performed by combining together 
ground-based radio telescopes with the orbiting 10-meter an¬ 
tenna of the SRT. The data are recorded at a rate of 128 Megabits 
per second, providing a total bandwidth of up to 32 MHz per cir¬ 
cular polarization channel ( Andreyanov et al.|2014[ l. The satel¬ 
lite telemetry and science data are received and recorded at track¬ 
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ing stations in Puschino ( [Andreyanov et aL]|2014| l and Green 
Bank ( [Ford et ar]|2014[ l and later transferred for correlation at 
one of the RadioAstron correlation facilities. Accurate recon¬ 
struction of the orbit and of the momentary state vector of the 
Spektr-R spacecraft is achieved by combining together radiomet¬ 
ric range and radial velocity measurements made at the satellite 
control stations in Bear Lakes and Ussurijsk (Russia), Doppler 
measurements obtained at each of the tracking stations, and laser 
ranging and optical astrometric measurements of the spacecraft 
position in the sky ( [Khartov et al.[2014{|Zakhvatkin et al.[2014[ l. 
The accuracy can be further enhanced by simultaneous VLBI 
measurements performed on the narrow-band downlink signal 
from the spacecraft ( [Duev et al.|2012| l. This was demonstrated 
( [Duev et al.|2015[ ) using the data from the ESP observations of 
0642+449 discussed here. 


The Spektr-R satellite has a highly elliptical orbit, with the 
apogee reaching up to 350,000 km and the perigee varying in the 
7,000-80,000 km range. The orbital period of the satellite varies 
between 8.3 and 9.0 days ( [Kardashev et al. [201 3[ [20145) 1. These 
factors make imaging observations with RadioAstron particu¬ 
larly challenging, with optimal conditions realized typically only 
for a limited fraction of sky and for limited periods of time. The 
imaging observations are planned by simulat ing the observing 
conditions in the FakeRaT software pac kag^^ ([Zhuravlev 2014 1 , 
which is based on the FakeSat software ([Murphy j 1991 [ [Murphy j 
let al.|1994K 


The Fourier spacings provided by the ground-space baselines 
to the SRT are typically limited to a very narrow range of posi¬ 
tion angle and often does not reach down to the spacings of the 
baselines between ground-based telescopes. In order to reduce 
the adverse effect of these two factors, imaging with RadioAstron 
is performed during the perigee passages of the SRT (perigee 
imaging), typically with 12-24 hour-long observing segments. 
Whenever feasible, these perigee imaging segments are aug¬ 
mented by several 1-2 hour-long segments (long-baseline track¬ 
ing) allocated during the following days and including a smaller 
number (3-4) of ground antennas. The long-baseline tracking 
provides information on largest Fourier spacings, enabling de¬ 
tecting the emission on the finest angular scales. 


2.1. Correlation of RadioAstron experiments 


Correlation of RadioAstron experiments is performed at three 
different correlator facilities; the ASC RadioAstron correlator 
( [Kardashev et al.|2013[ Andrianov et al.|2014[ l, the JIVE SFXC 
correlator ( Kettenis 2010 1, and the DiFX software correlator 
( [Deller et al.|2007| [201 l| l. In order to enable the correlation of 
RadioAstron data to be made in the DiFX correlator, the DiFX 
code has been upgraded by the Max-Planck-Institute for Radio 
Astronomy (MPIfR) in Bonn, addressing the specific aspects of 
data and telemetry formats of the SRT telescope on-board the 
Spektr-R spacecraft ( [Bruni[2014[|Bruni et al.[2015| l. 


The upgrade has comprised a new routine to read the native 
RadioAstron Data Format (RDF) into the DiFX correlator and 
a number of modifications to the correlator code enabling it to 
process data recorded at the SRT. The delay model server CALC 
from the CALC/SOLVE packag^used by the DiFX correlator is 
modified to be able to calculate delay information for telescopes 
with arbitrary coordinates and velocities, accounting for the gen¬ 
eral relativistic effects in the gravitational field of the Earth. The 


' http://www.asc.rssi.ru/radioastron/software/soft.html 
^ http://gemini.gsfc.nasa.gov/solve/ 



























































































Lobanov et al.: RadioAstron polarization image of 0642+449 


DiFX metadata system is upgraded to deal with the changing 
position and velocity of the spacecraft as a function of time. 

Satellite orientation parameters provided in the telemetry are 
used in the upgraded version of the code for calculating the 
equivalent of parallactic angle correction for the spaceborne an¬ 
tenna. The estimated accuracy of the orientation parameters is 
ss 1" ( |Lisakov et al.|2014l l. This enables accurate calibration of 
the instrumental polarization of the SRT. 

The timestamps for the SRT data (set by the tracking station 
clock at the moment it begins recording individual VLBI scans) 
are modified to refer to the relative time of arrival of the astro¬ 
nomical signal at the spacecraft. This is done by calculating the 
delay for the transmission of the signal from the spacecraft to the 
tracking station. These corrections are inserted to the metadata 
handler and to the delay model server, which provides the full 
compatibility of the RadioAstron data stream with data streams 
from ground antennas. 

The upgraded DiFX correlator is now used at the MPIfR to 
correlate RadioAstron data and being prepared for merging with 
the trunk version of the DiFX correlator code. 


3. RadioAstron observation of 0642-1-449 


The 1.6 GHz RadioAstron observation of 0642+449 (RadioAs¬ 
tron project code raesl2-, global VLBI project code GK047) 
was made on 9-10 March 2013, with three observing segments 
allocated for the experiment. The main perigee imaging seg¬ 
ment, GK047A (raesl2a; 10:00UT/09/03-01:00UT/10/03) was 
complemented by VLBI observations of the satellite down¬ 
link signal at 8.4 GHz (segment GK047B) made for the pur¬ 


pose of improving the Spektr-R orbit determination (Duev et al. 
[20T5] |. The third segment, GK047C (raesl2b; 14:00UT/10/03- 
15:00 UT/10/03) was introduced as a long-baseline tracking seg¬ 
ment. The SRT data and the satellite telemetry were recorded by 
the RadioAstron satellite tracking station in Puschino. Basic pa¬ 
rameters of the ground telescopes participating in the segments 
GK047A and GK047C are given in Table 

The data were recorded in two polarization channels (left and 
right circularly polarized, LCP and RCP), with a bandwidth of 
32 MHz per channel, centered at a frequency of 1.6601875 GHz. 
In each of the channels, the total bandwidth was split into two 
independent intermediate frequency (IF) bands of 16 MHz. 

The space VLBI segments on 0642+449 were observed in 
a 40/120 minute on/off cycle, with the gaps between the seg¬ 
ments required for cooling the motor drive of the on-board high- 
gain antenna of the Spektr-R. During these gaps, the strong ra¬ 
dio sources 1328+207 (3C 286) and 0851+202 (OJ 287) were 
also observed with the ground antennas as polarization calibra¬ 
tors. Fig.[2shows the resulting visibility coverage of the Fourier 
domain (uv coverage, expressed in units of spatial frequency de¬ 
fined as a ratio of instantaneous projected baseline length to the 
observing wavelength, d) of the data recorded for 0642+449. 


3.1. Data correlation 

The data correlation was performed at the DiFX software cor¬ 
relator of the MPIfR. Fringe search for RadioAstron was per¬ 
formed for every scan, and ad-hoc values for delay offsets and 
rates were determined and applied, in order to compensate for 
the acceleration terms of the spacecraft. The initial fringe-search 
was performed with 1024 spectral channels per IF and an inte¬ 
gration time of 0.1 sec (providing respective search ranges of 10 
microseconds and 10“® sec/sec for the delay and the delay rate). 



Fig. 1. Coverage of the Fourier domain (uv coverage) of the RadioAs¬ 
tron observation of 0642+449 at 1.6 GHz (d = 18 cm), plotted in units 
of Md. The central part corresponds to the perigee imaging segment 
(GK047A) and the short segment at v ~ 400 Md represents the long- 
baseline tracking segment (GK047C). Color marks the upper limit of 
brightness temperature obtained from visibility amplitudes (for details, 
see discussion in Sect.|4. l|l. 


in order to accommodate for potential large residual delays at 
the SRT. Fringes between the SRT and the reference antenna of 
Effelsberg were found with a residual delay offset of only about 
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Table 1. Radio telescopes participating in RadioAstron observation of 0642+449 


Telescope 

Code 

D 

[m] 

SEFD 

[Jy] 

Observing time [UT] 
GK047A GK047C 

^max 

[De] 

Spektr-R SRT (RU) 

RA 

10 

2840 

10:00-01:00 

14:00-15:00 


Effelsberg (DE) 

EF 

100 

19 

10:00-01:00 

14:00-15:00 

5.87 

Jodrell Bank (UK) 

JB 

76 

65 

10:00-01:00 

14:50-15:00 

5.93 

WSRT (NL) 

WB 

66f 

40 

13:00-01:00 

14:00-15:00 

5.90 

Torun (PL) 

TR 

32 

300 

10:00-01:00 


2.76 

Urumqi (CH) 

UR 

25 

300 

10:00-23:00 


3.02 

Shanghai (CH) 

SH 

65 

670 

10:00-19:00 


3.26 

Noto (I) 

NT 

32 

784 

10:30-01:00 


2.69 

Hartebeesthoek (SA) 

HH 

26 

430 

14:30-20:00 


1.91 

Green Bank (USA) 

GB 

100'' 

10 

17:00-01:00 


1.63 

Zelenchukskaya (RU) 

ZC 

32 

300 

19:30-23:30 


1.67 


Notes: D - antenna diameter (t - equivalent diameter); SEFD - system equivalent flux density (system noise), describing the 
antenna sensitivity; B^ax - largest projected baseline in the data, in units of the Earth diameter. 



Fig. 2. Total delay rates of the SRT obtained from baseline-based 
fringe fitting of the correlated data on the baseline to Effelsberg and ac¬ 
counting for the clock rates at each of the two telescopes. For each SRT 
observing segment, the rates are obtained from fringe fitting solutions 
with a time interval of 10 minutes. The rates are plotted separately for 
each IF and polarization channels. The plot reflects the accuracy of the 
velocity determination in the reconstructed orbit of the SRT. All of the 
velocity and acceleration values correspond to the baseline projections 
of the true residual offset, velocity, and acceleration of the spacecraft. 


1 microsecond or less, and a delay rate of about 2x10“" sec/sec 
(60 cm/s). The fringe solutions for the SRT were typically found 
at high signal-to-noise ratio (SNR): all but three scans showed 
fringes with SNR>10. For the three scans with low SNR, values 
of delay offset and rate were interpolated from adjacent scans. 
Final correlation was performed with 32 spectral channels per 
IF (with the corresponding channel width of 500 kHz) and 0.5 
sec of integration time. 


3.2. Post-correlation data reduction 


The correlated data were reduced in several steps using AIPS0 
and Difinap (Shepherd|1997 20111. The initial calibration was 


^ Astronomical Image Processing Software of the National Radio As¬ 
tronomy Observatory, USA 
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Fig. 3. Bandpass amplitudes and phases on ground-space baselines be¬ 
tween the SRT and Fffelsberg (top two rows) and Urumqi (bottom two 
rows), after applying global fringe fitting solutions. The bandpasses for 
total intensity (parallel hands: LL, RR) and polarized emission (cross¬ 
hands: LR, RL) are shown. 


performed in AIPS, imaging was done using Difmap, and both 
packages were used for calibrating the instrumental polarization 
and providing absolute calibration of the polarization vectors. 
The segments GK047A and GK047C were calibrated indepen¬ 
dently and combined together for the final imaging. Polarization 
calibration was performed using the segment GK()47A (the par¬ 
allactic angle coverage of GK047C was insufficient for using this 
segment for the polarization calibration). 
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Lobanov et al.: RadioAstron polarization image of 0642+449 



UV radius (10® A) 


Fig. 4. Visibility amplitude (top) and phase (bottom) distributions as a function of uv radius, overplotted with the CLEAN model (red color) 
obtained during the hybrid imaging of the source structure. 


The a priori amplitude calibration was applied using the de¬ 
fault antenna gains and the system temperature measurements 
made at each antenna during the observation. For the SRT, the 
sensitivity parameters measured during the IOC ( [Kovalev et al. 
|2014j l were used. The data were edited using the flagging in¬ 
formation from the station log data. Parallactic angle correction 
was applied, to correct for feed rotation with respect to the target 
sources. 


3.2.1. Fringe fitting 


The data were fringe fitted in two steps, first by applying manual 
phase-cal corrections and then using the global fringe search for 
antenna rates and single and multi-band delays. The group delay 
difference between the two polarization channels was corrected. 
The resulting post-fringe residual delay rates obtained for the 
SRT are plotted in Fig. These slowly evolving residuals agree 
well with the expected accuracy of the orbital velocity determi¬ 
nation for the SRT ( jKardashev et al.|2013 i, but they still should 
be viewed only as an indicator of the fidelity of the fringe fitted 
data, while a more detailed discussion of the accuracy of the or¬ 


bit determination for the SRT is presented elsewhere (Duev et al. 
[20T5] l. 


3.2.2. Bandpass calibration 

After the fringe fitting, the receiver bandpasses were corrected. 
The bandpass-calibrated data on the baselines between the SRT 


and the ground antennas in Effelsberg and Urumqi is shown in 
Fig. [3] The interferometric visibility signal is clearly detected 
both in the total intensity (LL and RR) and the cross-polarization 
(LR and RL) channels, demonstrating the excellent quality of the 
data. Table l^lists typical amplitude and phase errors of the vis¬ 
ibility data on the longest baselines to the SRT observed in the 
segment GK047A, demonstrating that the source was detected 
on baselines to all participating antennas (the errors were calcu¬ 
lated by vector averaging the calibrated data over two IFs and 
over a scan length of 10 minutes, and hence these errors are con¬ 
servative estimates containing also contributions from residual 
bandpasses and offsets between the IFs, and from amplitude and 
phase variations over the scan time). 


3.2.3. Polarization calibration 

After applying the bandpass calibration, the calibrated Stokes I 
data were exported to Difmap and imaged using hybrid imag¬ 
ing with CLEAN deconvolution and self-calibration (with full 
phase self-calibration and limited amplitude calibration using a 
single amplitude gain correction per antenna). The resulting self- 
calibrated visibility amplitudes and phases (and their CLEAN 
model representations) are shown in Eig. for all of the base¬ 
lines as a function of Eourier spacing {uv distance). The visibili¬ 
ties indicate that the source structure is clearly detected up to the 
longest ground-space baselines of the observation. 

The polarization calibration proceeded then with the self- 
calibrated model of the source structure obtained. The model 
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Table 2. Typical amplitude and phase errors on the baselines to the SRT 


Ant. 

Code 

Polarization 

LL,RR 

LR,RL 

^amp 

[%] 

n 

^ amp 
[%] 

^ph 

[°] 

EE 

5.5 

2.9 

22.6 

50.2 

JB 

9.3 

3.0 

30.7 

51.9 

WB 

6.0 

3.1 

30.5 

46.8 

TR 

11.6 

4.6 

38.3 

54.2 

UR 

13.6 

6.5 

45.0 

60.5 

SH 

53.9 

56.6 

53.9 

84.6 

NT 

18.3 

9.3 

54.1 

80.0 

HH 

52.8 

65.1 

50.8 

67.2 

GB 

6.5 

2.6 

33.4 

42.7 

ZC 

25.1 

13.8 

50.6 

57.7 


Notes: Columns represent average fractional amplitude, crj^p 
and phase, cTph, errors of the total (LL, RR) and cross-polarized 
(LR, RL) correlated flux density detected on the SRT baselines 
to the ground radio telescopes. 


Table 3. Instrumental polarization (D-terms) 


Antenna 

RCP 

LCP 

m 

[%] 

A 

[°] 

m 

[%] 

A 

n 

RA 

6.8+0.3 

82+6 

8.2+0.8 

98+1 


7.0+0.2 

81+8 

8.7+0.9 

98+2 

EE 

2.6+0.3 

10+2 

2.3+0.1 

167+9 


3.1+0.3 

8+3 

2.4+0.1 

185+9 

JB 

1.3+0.4 

140+7 

3.5+0.6 

336+3 


0.8+0.4 

306+7 

4.3+0.6 

330+3 

WB 

4.6+0.5 

176+3 

2.5+0.2 

166+13 


4.6+0.9 

182+3 

2.4+0.4 

199+12 

TR 

8.0+0.5 

160+3 

7.2+0.6 

14+8 


8.7+0.4 

164+3 

7.6+0.3 

31+8 

UR 

11.6+1.0 

271+12 

10.2+1.1 

56+18 


12.5+1.1 

294+12 

12.2+1.2 

93+18 

SH 

6.3+0.5 

173+31 

5.5+0.4 

13+30 


6.3+1.0 

170+32 

4.5+0.8 

349+30 

NT 

5.8+0.3 

247+5 

5.9+0.2 

239+6 


5.3+0.2 

184+5 

5.4+0.2 

178+8 

HH 

21.7+1.8 

151+23 

19.9+5.2 

80+35 


18.1 + 1.8 

198+23 

9.5+5.2 

149+35 

GB 

5.4+0.3 

248+6 

5.1+0.8 

295+10 


4.8+0.3 

236+6 

3.5+0.8 

315+10 

ZC 

8.0+0.7 

160+7 

5.2+1.0 

146+3 


9.5+0.7 

146+7 

7.2+1.0 

141+3 


Notes: For each antenna, listed are the fractional amplitude, m, 
and phase, x, of the instrumental polarization (polarization leak¬ 
age) in the right (RCP) and left (LCP) circular polarization chan¬ 
nel and in the first (IFl, top row) and second (IF2, bottom row) 
intermediate frequency channel. The D-term phases have been 
corrected for the phase offsets obtained from the EVPA calibra¬ 
tion described in Sect. 13.2.41 


was imported into AIPS and used as the calibration model for 
the full Stokes data, with amplitude and phase calibration ap¬ 
plied on time intervals of 1 minute. 
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Table 4. Effelsberg measurements of total and polarized flux densities 
of the targets 


Source 

S tot 

[Jy] 

S pol 

[Jy] 

A 

[deg] 

0642+449 

1.32 + 0.02 

0.022 + 0.008 

11.2 + 12.0 

0851+202 

2.77 + 0.03 

0.069 + 0.016 

32.1 + 1.7 

1328+307 

13.68 + 0.21 

1.37 + 0.10 

33.3+ 1.4 


Notes: Measurements were made on 9-10 March 2013. Column 
designation; 5 tot - total flux density; 5poi - polarized flux den¬ 
sity; m - percentage of polarization; x - position angle of the 
electric vector position angle (EVPA). 


The instrumental polarization (polarization leakage terms or 
D-terms, comprising the fractional amplitude and phase) was ob¬ 
tained using the LPCAL method ( |Leppanen et al.|[l995| l, with 
Effelsberg used as the reference antenna. The resulting polar¬ 
ization leakages determined through this procedure are listed in 
Table Errors of the D-terms were estimated from averaging 
the respective values determined for each of the calibrators and 
for the target itself To estimate the errors for the SRT (which 
did not observe the calibrators), the D-terms obtained from each 
of the calibrators with the ground antennas were applied to the 
data on 0642H-449 and respective solutions for the D-terms of 
the SRT were obtained and compared with the solution obtained 
directly from the target data (hence the D-terms errors for the 
SRT might be underestimated). The amplitude of instrumental 
polarization of the SRT is found to be within 9% and with ex¬ 
cellent phase consistency between the lEs. This figure cannot 
be directly compared with the results from pre-launch measure¬ 
ments, as those measurements were made only for the antenna 
feeds of the SRT, yielding an a; 2% leakage at 1.6 GHz ( |Tury- 
|gin| 2014) . The polarization leakage of the SRT is comparable 
to the instrumental polarization obtained for the ground anten¬ 
nas, demonstrating robust polarization performance of the SRT 
and ensuring reliable polarization measurements with data on 
all ground-space baselines. Similar estimates of the L-band po¬ 
larization leakage were obtained for the SRT from the analysis 
of the RadioAstron AGN survey observations (|Pashchenko et al.| 

[20l3] l. 

After applying the D-terms determined for the SRT and the 
ground antennas, the calibrated full Stokes data were exported 
into Difmap, imaged in Stokes I, Q, and U, and combined to¬ 
gether, providing images of the total intensity, linearly polarized 
intensity, and polarization vectors (electric vector position angle, 
EVPA). 


3.2.4. EVPA calibration 

The absolute EVPA calibration was performed using single-dish 
measurements of the total flux density polarization of the tar¬ 
get and calibrator sources. These measurements (listed in Ta¬ 
ble 1^ were made concurrently with the RadioAstron observa¬ 
tion at the Effelsberg 100-meter telescope. The averaged EVPA 
of the polarized emission detected in the VLBI data was related 
to the polarization vectors that were measured with the Effels¬ 
berg antenna, and the resulting EVPA offsets were applied to the 
VLBI images. The resulting ground array images of the calibra¬ 
tor sources are shown in Eig. and the image parameters are 
listed in Table |5] 





































Lobanov et al.: RadioAstron polarization image of 0642+449 


Linearly Polarized Intensity (mJyAteam) 

1.09 4.71 8.32 11.94 15.56 



6 4 2 0 -2 

Relative Right Ascension (mas) 


Fig. 6. RadioAstron space VLSI image of 0642+449 of total (contours) and linearly polarized (color shades) emission, made using the uniform 
data weighting. The orientation of linear polarization is indicated by vectors. Basic parameters of the total intensity image are listed in Table 
Labels indicate the modelfit decomposition of the source structure described in Table|^ 


4. Polarization image of 0642+449 

The RadioAstron polarization image of 0642+449 is shown in 
Fig.and the basic parameters of the image are listed in Ta¬ 
ble The source shows a compact core-dominated structure, 



Fig. 5. Ground array images of the calibrator sources 1328+307 (left) 
and 0851+202 (right), made using the uniform data weighting. Im¬ 
age parameters are listed in Tahle|^ Contours show the total intensity; 
colors indicate linearly polarized intensity; vectors show the EVPA ori¬ 
entation. 


with a resolved nuclear region which contains about 98% of the 
total flux density in the image. The weak and extended jet fea¬ 
ture, Jl, located at about 4 mas separation from the nucleus has 
a peak brightness of a; 15mJy/beam, which is not surprising 
considering energy losses due to adiabatic expansion of the flow 
up to the estimated deprojected linear separation of 2.3 kpc of 
this feature from the core of the source (assuming the jet view¬ 
ing angle of 0.8°; Pushkarev et al. |2009|). The linear separation 


can be a factor of ~ 3.5 smaller, if the observed misalignment 
of ss 2° between the position angles of C2 and Jl reflects the 
physical change of the jet direction. The observed concentra¬ 
tion of the polarized flux density along the southern edge of this 
feature may reflect interaction with the ambient medium occur¬ 
ring in this region, although it may also result from the errors of 
the D-terms determination or from other instrumental effects (cf. 
[Gabuzda & G6mez|2001[rHovatta et al.|2012| l. 

The core region is clearly resolved and extended along a P. A. 
of ^ 81°, It can be decomposed into two circular Gaussian com¬ 
ponents. The resulting three-component decomposition of the 


Table 5. Parameters of total intensity images 


Source 

S tot 

“^peak 

S neg 

^rms 

Beam 

0851+202 

2.06 

1.73 

-19.0 

6.5 

12.2,2.93,12.6 

1328+307 

7.87 

0.89 

-47.1 

8.6 

7.48,2.73,14.9 

0642+449 

1.21 

0.81 

3.6 

0.7 

1.38,0.81,66.5 


Notes: Column designation; 5tot[Jy] - total flux density; 
5peak [Jy/beam] - peak flux density; Sneg [mJy/beam] - maxi¬ 
mum negative flux density in the image; cr^jns [mJy/beam] - rms 
noise in the image; Beam: major axis, minor axis, position angle 
of major axis [mas, mas, °]. 


Article number, page 7 of 


10 

























































Table 6. Decomposition of compact structure in 0642+449 


Comp. 

‘^tot 

[mJy] 

r 

[mas] 

n 

e 

[mas] 

^min 

[mas] 

Tb 

[10" K] 

^ pol 

[mJy] 

^^pol 

[%] 

Xpol 

n 

Cl 

462+16 

0.38+0.01 

-88.3+0.7 

0.40+0.01 

0.16 

12.8+0.6 

7.7+0.8 

1.7+0.2 

79+3 

C2 

721+19 

0.32+0.01 

70.1 + 1.1 

0.60+0.01 

0.13 

8.8+0.3 

9.9+0.7 

1.4+0.1 

89+2 

Jl 

20+9 

3.80+0.20 

92.0+3.2 

0.89+0.43 

0.75 

0.11+0.04 

1.7+0.7 

8.5+5.2 

39+11 


Notes: Gaussian model description; 5tot - total flux density; component position in polar coordinates (r, 0) with respect to the map 
center; component size, 0; minimum resolvable size, 0mm, of the respective component; brightness temperature, Tb, in the observer’s 
frame, derived for the parameters of the Gaussian fit. Polarization properties (measured from the Stokes U and Q images: 5poi - flux 
density of linearly polarized emission; mpd - fractional polarization; Xpoi - polarization position angle in the sky plane (measured 
North-through-East). 


source structure is presented in Table The two nuclear com¬ 
ponents, Cl and C2, are separated by 0.76mas (a; 410pc, de- 
projected). Similar structural composition of the nuclear region 
is also detected in the ground VLBI from the MOJAVE survey 
at 15 GHz ( Lister et al.|2013 1. The same structure is also present 
in the images of 0642+449 made at 24 and 43 GHz (|Fey et al.| 
[ 2000 ) 1 . 

Overall, the structure observed in 0642+449 may be similar 
to the main features of the jet in M 87 (cf. |Sparks et al.|1996| l, in 
which a prominent standing shock (component HST-1) is located 
at 200-280 pc (considering the uncertainty of the jet viewing an¬ 
gle) and a prominent feature A located at 2.8-3.7kpc distance 
from the core is believed either to be an oblique shock ( [Bicknell] 
& Begelman|1996]l or to result from the helical mode of Kelvin- 


Helmholtz instability ( [Lobanov et al.|200^ . This may also ex¬ 
plain the apparent displacement between the total intensity and 
linearly polarized emission in the region Jl, which is similar to 
the polarization asymmetries observed in M87 at these linear 
scales ( jOwen et al.|1989| l. 

The polarized emission detected in 0642+449 is dominated 
by the two nuclear components, with fractional polarization 
reaching 15.6 mJy (^ 2%), and the magnetic held showing pre¬ 
dominantly transverse orientation with respect to the general 
jet direction (assuming that the innermost component Cl cor¬ 
responds to the canonical “VLBI core” region in which the 
emission becomes optically thin at the frequency of observa¬ 
tion; |^nigJT98T|l. The polarization vectors vary by about 10° 
across the nucleus, with the polarization angle observed in the 
most compact region Cl (=s 80°) coming close to the polariza¬ 
tion angle measured with Effelsberg. A substantial rotation of 
the EVPA vectors is seen in some of the MOJAVE images of 
0642+449, which is likely to result from plasma propagation in 
a jet observed at an extremely small viewing angle. The weak 
(~ 2mJy/beam) linearly polarized emission in the jet compo¬ 
nent Jl corresponds to a fractional polarization of ~ 15%, with a 
large error margin on this estimate. The magnetic held is roughly 
longitudinal in this region. 

The MOJAVE databas^ lists, for the frequency of 15 GHz, 
an EVPA of 122° in March 2011. Analysis of the multifre¬ 
quency MOJAVE observations ( Hovatta et al.|201^ yields, for 
0642+449, a median rotation measure of -268rad/m^ over the 
whole source structure, with a slightly higher rotation measure 
of -280 rad/m^ over the compact core region. 

Assuming a linear trend in the EVPA in time as derived 
from the two latest MOJAVE observations, the polarization an¬ 
gle should be ~ 105° in March 2013 at 15 GHz. We measure 
X - 80°-90° over the core region. Eormally, this corresponds to 


http://www.physics.purdue.edu/MOJAVE/ 
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a rotation measure of -15 rad/m^, under assumption that the ro¬ 
tation measure is not affected by the core shift. It would then re¬ 
quire about three turns of phase to get this value to be reconciled 
with the RM reported from the MOJAVE measurements. Alter¬ 
natively, the smaller RM estimated from the RadioAstron mea¬ 
surements may reflect the reported trend of decreasing RM with 
decreasing frequency ( Kravchenko et al.|2015| l. It would still re¬ 
quire about one turn of phase to bring the estimated RM value to 
a better agreement with the empirical dependence |RM|(v) oc y* ® 
derived injKravchenko et ^(|2015]l. 


4.1. Brightness temperature 


The Gaussian modelflt described in Table |6] can be used for 
obtaining estimates of brightness temperature (also presented 
in the same Table), with the highest brightness temperature of 
1.3 X lO'^K associated with the most compact core feature Cl. 
This value is somewhat lower than the one that should be ex¬ 
pected at 1.6GHz, based on the trend log[Q7’b[K] = 13.3 - 
1.1 logjQy[GHz] that can be obtained from the existing mea¬ 
surements at 5, 15, and 86 GHz (Lobanov et al.|2000) [Kovalev) 


jet al. [2005] [Dodson et al.|2008HLee et al.|2008] l. However, since 

the morphology of the emitting region can be different from a 
simple Gaussian shape, the actual brightness temperature of the 
emission may differ from the modelflt-based estimates, for in¬ 
stance, if the brightness distribution is smooth and the respective 
brightness temperature is largely determined mostly by the trans¬ 
verse dimension of the flow ( [Lobanov[2015[ l. To investigate this 
possibility, the modelflt-based estimates of Tb can be compared 
with estimates obtained directly from visibility amplitudes and 
their errors, which yields the absolute minimum brightness tem¬ 
perature, Tb^min, and an estimate of the limiting brightness tem¬ 
perature, Tb.Mm that can obtained from the data under the require¬ 
ments that the structural detail sampled by the given visibility is 
resolved ([Lobanov[2015[l. 


These two estimates are compared in Eig.jTjwith the bright¬ 
ness temperatures calculated from the Gaussian modelflt. This 
comparison shows that the visibility amplitudes on the space 
baselines longer that « 140 Md require the brightness temper¬ 
ature to be larger than 10'^ K, which agrees with T\, obtained 
for the most compact feature Cl of the Gaussian modelflt. The 
limiting brightness temperature = 3.6 x 10'^ K can be es¬ 
timated at the longest space baselines. These visibilities con¬ 
centrate however only within a narrow range of position angles 
(-10/2 + 0/1), with corresponding structural sensitivity limited 
to P.A. ^ 80°, which essentially coincides with direction of the 
inner jet (hence resulting in lower visibility amplitudes and, con¬ 
sequently, lower estimates of Tbjim)- 











































Lobanov et al.: RadioAstron polarization image of 0642+449 



Fig. 7. Visibility-based estimates of the brightness temperature in 0642+449, calculated in the observer’s frame. Data points represent the 
minimum (dark gray), rb,min> and maximum limiting (light gray), Lbun,, brightness temperatures derived from individual visibilities measured 
at different uv distances, q jLobanov|2015| l. The two-dimensional distribution of rbji,„ is also presented in Fig. [T] with a color wedge. Braces 
above and below the brightness temperature data indicate the range of uv distances covered by the space baselines ot the RadioAstron observation. 
The histograms show the respective brightness temperatures averaged over bins of 5 MT in size. The thick dashed line indicates the brightness 
temperature, Tb^mod = 1-3 x 10*^ K, derived for the most compact feature. Cl, of the Gaussian fit of the source structure described in Table 
The dotted lines indicate the respective brightness temperatures of the other two Gaussian components. The minimum brightness temperature, 
Tb.min = 1-2 X 10'^ K, is constrained by the visibilities with q > 380 MT. The maximum limiting brightness temperature Tb nn, = 2.4 x 10*'* K can 
be estimated from the data on the SRT baselines to Noto and Hartebeesthoek, although potential problems with amplitude calibration cannot be 
ruled out. A more conservative estimate of Tbaim = 3.6 x 10*^ K can be obtained from the data on the longest space-ground baselines. 


Indeed, some of the Tb.Km estimated from the visibility am¬ 
plitudes at shorter space baselines and different position angles 
reach above lO'^ K, notably with a large spread of the estimates 
obtained at similar spatial frequencies. The observed spread re¬ 
flects the difference of signal-to-noise ratios (SNR) in the mea¬ 
surements on different baselines, and it also may indicate po¬ 
tential biases due to a priori gain calibration and amplitude self¬ 
calibration. The three particularly outstanding clusters of points 
in Fig.0 which have Tb^im > 10'"^K, come from visibilities on 
the basdines to Westerbork, Torun, Noto, and Hartebeesthoek. 

As the visibility SNR on these baselines is typically lower 
than on baselines to Effelsberg or Jodrell Bank (which do not 
yield such high Tb^iim at similar radial distances and position 
angles), it is most likely that the exceptionally high brightness 
temperature limits are caused by problems with the amplitude 
calibration of the antenna gains or specihc observing scans with 
these antennas. 

If the outliers are ignored and only the bulk of the Tb^iim es¬ 
timates are considered in the range of 100-200 MT, the average 
Tb.iim increases to 9.1 x 10*^ K, and this can be viewed as a viable 


upper limit on the brightness temperature of the most compact 
structure in the jet of 0642+449. This corresponds to brightness 
temperature of 4.0 x 10^^ K in the rest frame of the source. If the 
brightness temperature is generally determined by the transverse 
dimension of the flow, as suggested by the analysis of the MO¬ 
JAVE visibility data ( |Lobanov|2015| l, the inner jet of 0642+449 
should have a width of ^ 0.15 mas to satisfy the limiting bright¬ 
ness temperature estimated above. 


5. Summary 

The RadioAstron Early Science observations of the high-redshift 
quasar TXS 0642+449 have tested the polarization performance 
of the space radio telescope (SRT) at 1.6 GHz (A - 18 cm) 
and provided a framework for establishing a set of procedures 
for correlation, post-processing and calibration of space VLBI 
polarization data. The instrumental polarization of the SRT 
is found to be less than 9%, which enables robust reconstruc¬ 
tion of the polarization signal on the ground-space baselines of 
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RadioAstron observations and provides a basis for high-fidelity 
imaging of polarized emission. 

The total and linearly polarized emission from the quasar 
0642H-449 is detected on the SRT baselines of up to 75,560 km 
(5.93 Earth diameters) in length, corresponding to angular scales 
of 0.5 milliarcseconds (mas). The resulting RadioAstron image 
of linearly polarized emission in 0642H-449 has a resolution of 
0.8 mas, which is ~ 4 times better than the resolution of ground 
VLBI images at 18 cm. The image probes the physical condi¬ 
tions in the jet on angular scales as small as « 0.2 mas (linear 
scale of ss 1.5 pc) taking into account the SNR-driven resolution 
limits of the observation. 

The structure of total intensity and linearly polarized emis¬ 
sion recovered from the RadioAstron image of 0642H-449 sug¬ 
gests that this object is likely to be a “canonical” quasar with 
a powerful relativistic jet observed at an extremely small angle 
of sight and featuring a bright “core” of the jet (component Cl) 
and a likely recollimation shock (component C2) located at a 
deprojected distance of about 400pc away from the core. Each 
of these regions is only weakly polarized, with fractional polar¬ 
ization not exceeding « 2% and polarization vector implying a 
predominantly transverse magnetic field (in agreement with ex¬ 
pectations for strong shocks in jet). The weak and more distant 
feature J1 observed at a separation of 3.8 mas (=s 2kpc, de¬ 
projected) has a much stronger polarization, possibly exceeding 
a 10% level. The dominant magnetic field components is likely 
to be poloidal in this region, implying that it either represents a 
weak shock or may even result from plasma instability develop¬ 
ing in the jet. 

The RadioAstron data have also been used for making bright¬ 
ness temperature measurements based on the Gaussian decom¬ 
position (modelfit) of the structure and on estimates made di¬ 
rectly from the visibility data. The most robust estimate pro¬ 
vided by these measurements yields a brightness temperature of 
ss 3.6 X 10'^ K for the most compact region in the jet. Eurther 
analysis of the RadioAstron data indicate that the the brightness 
temperature of the radio emission from this region cannot be 
lower that ~ 10'^ K and is not likely to exceed ^ 9 x lO'^K 
(corresponding to=s 4xlO'^Kin the reference frame of the 
source). 
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